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In this work, Rydberg energy levels of the triatomic hydrogen molecule (H3) are studied with mul- 
tichannel quantum-defect theory. We extract the body-frame p-wave quantum defects from highly 
accurate ab initio electronic potential surfaces and calculate the quantum defects of higher angular 
momentum states in a long-range multipole potential model. Laboratory-frame quantum defect 
matrices emerge from a rovibrational-frame transformation carried out with accurate rovibrational 
states of Hj. Finally, we use the laboratory- frame quantum defects to calculate Rydberg energy 
levels for the fundamental neutral triatomic molecule H3. 
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I. INTRODUCTION 

The triatomic hydrogen molecule (H3) plays an impor- 
tant role in astrophysics because its cation Hj acts as a 
proton donor in chemical reactions occurring in interstel- 
lar clouds. As the simplest triatomic neutral molecule, 
H3 also attracts fundamental interest. Ever since its 
emission spectra were first observed by G. Herzberg in 
the 1980s [l|, H 3 has been studied extensively. Herzberg 
and co-workers measured infrared and visible emission 
spectra of H 3 in discharges through hydrogen and as- 
signed them to Rydberg-Rydberg transitions between 
n = 2 and n = 3 electronic states using empirical fits 
Helm and co-workers investigated the higher Ryd- 
berg states and ionization potentials of H3 by analyzing 
the photoabsorption spectrum [2j. In 2003, building on 
previous work of Schneider, Orel and Suzor-Weiner [3j, it 
was shown 0, [j| that intermediate Rydberg states of H3 
play an important role in the dissociative recombination 
(DR) process, H+ + e" -> H3 -> H 2 + H or H + H + H. 
Prior to the study of Ref. [j, [f| , the large discrepancy 
between the DR rate determined by experiment and pre- 
vious theory had not been resolved. Refs. [3, H| found 
that Jahn- Teller effects in H3 neglected in previous the- 
oretical studies couple the electronic and nuclear degrees 
of freedom and generate a relatively high DR rate via 
intermediate p-wave Rydberg-statc pathways. A recent 
alternative formulation developed by Jungcn and Pratt 
provides supporting evidence for this interpretation Q. 
Vervloet and Watson improved both the experimental 
techniques and empirical fits and reinvestigated the low 
Rydberg states that G. Herzberg had observed Q. Here 
we undertake an analysis of the Rydberg states with ab 
initio theory. One of the most successful techniques in 
treating Rydberg states by ab initio theory is multichan- 
nel quantum defect theory (MQDT) [1,0. Earlier studies 
[3, []jM3 have utilized MQDT to successfully describe 
the DR process. 

The application of MQDT to study molecular Rydberg 
energy levels treats the H 3 molecule as a Rydberg elec- 
tron attached to a Hjj" ion. The interaction between the 
Rydberg electron and the ion core is described through 
a smooth reaction matrix K or quantum defect matrix 
fi. K and /1 are simply related, e.g., for a single-channel, 



K = tan (71-^). We extract a body- frame reaction ma- 
trix from ab initio electronic potential surfaces for p- 
wave Rydberg states and calculate the body-frame reac- 
tion matrix for Higher angular momentum states (I > 1) 
by using the long-range multipole potential model. For 
the higher angular-momentum states, we neglect short- 
range interactions due to the nonpenetrating nature of 
the high / states. Here I denotes the quantum number 
of the Rydberg electron orbital angular momentum. We 
then construct the total laboratory-frame reaction matrix 
K through a rovibrational frame transformation, obtain- 
ing 



K w = 



(i\ a)K aa , (0/ 



(1) 



Here, Kai is the laboratory-frame reaction matrix ele- 
ment between the laboratory- frame eigenchannels \i) and 
\i'), and K aa > is the body-frame reaction matrix ele- 
ment between body-frame eigenchannels \a) and \a'). 
The rovibrational frame transformation is specified by 
the unitary transformation Ui a = (i \a). 

The process of constructing the rovibrational transfor- 
mation is similar to that described in Ref. 0and is based 
on the rovibrational wave functions of H^ . To calculate 
them, there are two important approximations adopted 
in Ref. [||, the rigid rotator approximation and the adi- 
abatic hypcrspherical approximation. The nonadiabatic 
coupling between different adiabatic hyperspherical chan- 
nels was included in later studies by using the slow vari- 
able discretization (SVD) approach in Ref. [H, In 
the present study, we abandon the rigid rotator approx- 
imation and consider the Coriolis interaction. In this 
way, we obtain very accurate rovibrational energy levels 
and wave functions of H^ that allow us to construct the 
rovibrational transformation. 

After the rovibrational transformation described by 
Eq.([T]) is carried out, we obtain the laboratory-frame 
K matrix and calculate the eigenenergies E of the H3 
molecule by solving the secular equation flo| 



det I tan (wu) + K\ = 0, 



(2) 



whe 



re v is a diagonal matrix with elements va = 
Here ET V denotes the ith rovibrational 



i/V2T^ 

energy level 



E). 
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The remaining of the article is organized as follows. 
Section 2 describes the detailed calculation of rovibra- 
tional states of H^" and shows how to use them to con- 
struct the rovibrational transformation. Section 3 de- 
scribes the calculation of the p-wave energy levels of H3 
using ab initio quantum defects. Section 4 discusses the 
long-range multipole potential model for higher angular 
momentum Rydberg states, and Section 5 gives our con- 
clusions. 



II. ROVIBRATIONAL STATES OF H+ 

In the cation Hjj~, three protons interact with each 
other under a potential surface. The potential surface 
was created by Refs. and is sub-micro-hartree accu- 
rate. With this potential surface, we solve the three-body 
Schrodingcr equation in the adiabatic hyperspherical rep- 
resentation. We use the SVD method to include the 
nonadiabatic coupling between different adiabatic chan- 
nels. We also include the Coriolis interaction that couples 
rotational angular momentum with vibrational angular 
momentum. 

After separation of the center-of-mass motion, the 
three-body system can be described by six coordinates, 
the Euler angles a, /3, and 7 and the three interparticlc 
distances ri2, ^23, and r$i- Alternatively, the interpar- 
ticle distances can be expressed in terms of hyper- 
spherical coordinates, i.e., the hyperspherical radius R 
with hyperspherical angles 9 and tp [13, HH , where 



n 2 = 3~ 1/4 P[1 + sin6>sin (p - tt/6)] 



1/2 



r 23 = 3~ 1/4 P [1 + sin 9 sin (tp - 5vr/6)] 



1/2 



+ sin 9 sin (<p + it/2)] 



1/2 



(3a) 



(3b) 



(3c) 



The hypcrangle 9 ranges from to 7r/2, while the hy- 
perangle tp ranges from to 27r. The hyperradius R ex- 
tends from to infinity. In hyperspherical coordinates, 
the Schddinger equation is given by 



1 



A 



2fiR 5 dR dR 2^K 2 



V(R,6,tp) 



(4) 

The Coriolis interaction is included in Eq. Q through the 
squared "grand angular-momentum operator" A 2 which 
is given by 



2^R 2 



To + T v c + T r , 



where, 



Tn 



2 9 ■ 99 9 
^R 2 sin 29 89 ^ W 



(5) 



(6) 



T 



<pC 



d 



/iP 2 sin 9 V dp 



cos ■ 



and 
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■J 



nR 2 (l-sm9) ,iR 2 (1 + sin(9) 2/iP 2 ' 



(7) 



(8) 



The operators (J x , J y , J z ) are the body- frame compo- 
nents of the total angular momentum of the ion, and 
/j = m/-\/3 is the three-body reduced mass of the sys- 
tem, where m is the atomic hydrogen mass. Solving the 
Schddinger equation Eq.Q directly should in principle 
provide accurate rovibrational energy levels, but it would 
require extensive computational time and memory to di- 
agonalize the full Hamiltonian matrix. Instead, we break 
the problem into two steps: first solve the hyperangular 
Schddinger equation in the adiabatic representation, and 
then later include the nonadiabatic coupling through the 
SVD method. 



A. Adiabatic representation 

In the adiabatic representation, the adiabatic poten- 
tials and channel functions are defined as solutions of 
the adiabatic eigenvalue equations, 

A 2 1 c 1 

V(R,6,tp) 



2\iR 2 S^R 2 



S„ (0; R) = U v OR) (O; R) 
(9) 



whose solutions depend parametrically on R. For each 
P, the set of $„ (f2; R) is orthogonal, complete, and can 
be used as a basis to expand the whole wave function, 
eg-, 



% (R, Q) = J2 R^' 2 K (R) ®u (fi; P). 



(10) 



Here we use f2 to denote the Euler angles and the two 
hyperspherical angles. To solve Eq.® numerically, we 



expand <!>„ (fi; R) in a set of basis <j> im m K + I , such that 



N+m+g, 



jm 2 I 

-N + m+ gi 



(11) 



jm,2 



where j, mi, N + , K + and gi are quantum numbers 
labeling basis functions in different degrees of freedom. 



^ I jm^K+ I sa ti snes the permutation symmetry of a three- 
fermion system, e.g., 



-^N + m + g I _ 



and 



A<f> 



-N + m+gi 
~jm 2 K+ ' 



where 

A = 1 — Pi 2 — P23 — P31 



P12P31 + P12P 



23- 



(12a) 
(12b) 
(13) 



The explicit form of the basis functions ^ m ^+' is given 
in Appendix A. 
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B. Slow variable discretization 

With the help of the discrete variable representation 
(DVR) basis ir n (R) , effects of nonadiabatic couplings be- 
tween different adiabatic potentials are included. The 
DVR basis has the following properties: 



dRir n (R) U v (R) tt„> (R) « U v (R n ) S n 



(14) 



where R n are the quadrature abscissas used to generate 
the DVR basis. The SVD approximation expresses the 
expansion of Eq. (fTU)) as 



^ {R, fl)=J2 c^R~ 5 / 2 TT n (R) (Q; R n ), (15) 



and rewrites Eq.@ as 



7 j x nn f yy 



+ [U v (R n ) - E™] e nv = 0, (16) 



where 



dRlT n {R) 



2 1 



1 d 
%l~dR2 



ir„, (R) dR, (17) 



and O nv ^ n i M is the overlap matrix given by 

o nu , n .„ = (^ v (n-,Rn)\ * M (n ; <)) 



(18) 

Finally, Eq. ((T6|) is solved for the expansion coefficients 
c l nv and E" '. The total rovibrational wave function is 
therefore given by 

nv jm 2 K+ 

(19) 

corresponding to the eigenenergy E™. Here i is the set 
of good quantum numbers N + m + v + giH + , where v + de- 
notes the vibrational quantum numbers. 



C. Accuracy of rovibrational energies of 

Next we compare our theoretical rovibrational en- 
ergy levels E™ of H^ with experimental energy lev- 
els Adopting the notation used in Ref. fl9| . we 
label the rovibrational states i by quantum numbers 
(N + , G){vi, v l 2 2 }(l\u). vi is the symmetric-stretch vibra- 
tional quantum number, v 2 denotes the quantum number 
of the asymmetric-stretch mode, l 2 describes the quan- 
tum number of the vibrational angular momentum, and 
G = \K + — Z2 1 - The fact that G instead of K + is a good 
quantum number implies that the Coriolis interaction 
couples rotational and vibrational angular momenta and 
makes levels with the same G nearly degenerate. How- 
ever, for levels with l 2 + and (N+ - \l 2 \) > G > 1, 
the degeneracy breaks, and we u (or I) to denote the up- 
per (or lower) energy level; these levels with a u or an I 
cannot be described by rigid rotator approximations. 



TABLE I: Comparison of several calculated rovibrational en- 
ergy levels of with experimental results [l9| . Only states 
with N + < 3 for the {0, 0°} and {0, l 1 } bands are shown here. 



Q.N. 


&cal 


rp c 


Q.JN. 




771 C 
J-^exp 




(cm ) 


(cm ) 




(cm ) 


(cm ) 


(l,l){0,0 u } 


64.128 


64.121 


(2,3){0,1 1 } 


2614.034 


2614.270 


(1,0){0,0 U } 


86.960 


86.960 


/r» o\ fn ill 

(2,2){0,1 } 


2723.708 


2723.962 






1 fSQ 9Q c i 




97^ ^1 ^ 


97^ 5R5 


(2,1){0,0 } 


237.335 


237.356 


(2,l){0,l 1 }u 


2790.086 


2790.344 


(3,3){0,0 } 


315.317 


315.349 


(2,0){0,1 1 } 


2812.504 


2812.850 


(3,2){0,0 } 


427.974 


428.018 


(S^HO.l 1 } 


2876.566 


2876.847 


(3,1){0,0 } 


494.712 


494.775 


^HO,! 1 }/ 


2931.091 


2931.366 


(3,0){0,0 } 


516.823 


516.873 


(3,2){0,1> 


2992.151 


2992.436 


((MHO,! 1 } 


2521.183 


2521.411 


(3, 1){0, l 1 }/ 


3002.348 


3002.905 


(MHO,! 1 } 


2547.996 


2548.164 


(3,0){0,1 1 } 


3025.663 


3025.941 


(MHo,! 1 } 


2609.302 


2609.541 


(S,!)^,! 1 }^ 


3063.181 


3063.951 



a The quantum numbers labeling the energy levels, in the 
notation (N + , G) {vi, w 2 2 }(Z|«) described in the text. 
b Theoretically calculated results from this work. 
Experimentally determined energies from Ref. 



Table 2 compares the rovibrational energy levels, 
calculated for N+ < 3 states of {0,0°} and {0, l 1 } 
bands, with the experimental results of Lindsay and 
McCall. [lj| The agreement is good, with a rms difference 
of 0.281 cm -1 for the levels shown in the table. Higher- 
rovibrational energy-level calculations also exhibit good 
agreement. For our calculated energy levels up to around 
9000 cm -1 with N + < 4, the rms difference between our 
calculation and the experimental results of Ref. [l9| is 
0.657 cm -1 . 



D. Rovibrational- frame transformation 

Next we describe in detail how to construct the 
rovibrational-frame transformation using the ionic rovi- 
brational eigenstates. In the laboratory frame, the 
H$ +e~ system is described by the electron orbital angu- 
lar momentum I and its projection A onto the laboratory 
z-axis, and by A + , m + , gj, and parity of the ion core. 
Hence, we construct the wave function of the + e~ 
system as a sum of products of the ionic rovibrational 
wave function and electronic wave function, of the form 



^irS (R,n)Y lx (e e ,cp e ), 



(20) 



where 9 e and tp e are spherical angles of the electron in 
the laboratory frame. The radial part of electronic wave 
function is not shown here since we apply the MQDT 
method for that degrees of freedom. Specifically, we want 
to construct the laboratory eigcnchanncl function \i) with 
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a definite laboratory-frame total angular momentum N: 
\i) = | N + , v +^{ N > l > m > n <9-0_ j n ^he position representa- 
tion, it takes the form as, 



N+ ,m+;L\ ' N+m+ / ' A ' 



Y l} 



(21) 



with some manipulation of algebra, where 9' e and ip' e are 
the spherical angles of the electron in the body-frame. 
Using Eq. (|24|) . we derive the transformation between 
the body- frame and laboratory-frame states as, 



where C N + l +.> \ denotes the appropriate Clebsch- 
Gordan coefficient. 

In the body-frame, the is specified state by the pro- 
jection of theelectron angular momentum on the molec- 
ular Z-axis A, by the total angular momentum of the 
system N including the electron contribution I and the 
projection m of iV on the laboratory z-axis. Applying 
the transformation between body-frame and lab-frame 

Y n {6 e ,<p e ) = J2 [D[a (a, P, 7)] * *ia (22) 



the expansion of the product of two Wigner functions 

~,N,K s~<N,m 

t 1J mK { - 

N 

we have the following equation, 



U m+K+ U \A 



U mK { - J N+ ,K+\l.A^ N+m+;l,\> 



(23) 



(a \i) = (R,Q;A \ N+ , u+) (»^) = 
where ^^+^+a is defined as, 



(25) 



N+v+A 



E 



ct + K. t (R) 



X E «SU+ ( R ^-Z% +K+ A^) 



\Nmgi 



A 

E(- 



N,rr 
N+, 



:I.X 



n N 

U m+K+ 



Y lx (6 e ,(p e ) (24) 



;-a r N + ,K+ 
°i,-A;Ar,if 



[D? nK } YiA^cp'J, 



The explicit form of ^^n+k+a 1S §i vcn m Appendix 
A. The rovibrational frame transformation can be ac- 
complished as follows: 



K 



(Nm gi n+) _ 
N+v+;N+'v + ' 



f dQdn E y 

A, A' 



'*K A a> (Q)^ NmII+91 



N+v+A 



N+'v+'A" 



(27) 



r 



where Q,e denotes the Euler angles, and Q contains the 
three vibrational degrees of freedom. 



where U denotes the unitary transformation that diago- 
nalizcs the laboratory-frame quantum-defect matrix. 



III. P-WAVE ENERGY LEVELS OF H 3 



A. Body-frame quantum defects for p-waves 



The body-frame reaction matrix for a p-wave elec- 
tron is described by the short-range interaction extracted 
from an ab initio calculation |2(|. In practice, the quan- 
tum defects are smoother than the reaction matrix el- 
ements because the latter can have poles. Hence in 
this work, we extract the body-frame quantum defects 
M (Q) from the ab initio energy surface directly. After 
replacing A'aa' (Q) by /xaa' (Q) in Eq. (|2~7} . we perform a 
rovibrational transformation to get the laboratory-frame 
quantum-defect matrix. Finally, we get the laboratory K 
matrix by using the eigenvalues /i e from the laboratory- 
frame quantum-defect matrix, 



K = Utan(TTH e )U q 



(28) 



Because of Jahn- Teller effects, the body-frame 
quantum-defect matrix is generally not diagonal in the 
electronic projections A, A'. Similar to the body- frame 
K matrix proposed by Staib and Domcke [2l| , it has the 
form 



m(Q) 



Moo (C) 

Mn(Q) mi-i(Q) 
m-ii(Q) m-i-i(C) 



(29) 



We express each matrix element by using the vibrational 
symmetry coordinates Q = (Qi, Q x , Qy) as 

Qx = f^= (An + Ar 2 + Ar 3 ) , (30a) 
V3 
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Q x = f-±= (2Ar 3 - Ar 2 - An) , (30b) 

and 

Q v = f (An - Ar 2 ) . (30c) 

where / = 2.639 255 bohr -1 is a constant and Ar; de- 
scribe displacements of the nuclei from the equilibrium 
geometry at which r%2 = r 2 3 = r^i = r equ i = 1.6504 a.u.. 
For example, An = ?" 2 3 — ?"equi- (Qx,Qy) can be alter- 
natively described by another pair of coordinates (p, (f) 
as Q x = pcos4> and Q y = psin0. Q\ describes the sym- 
metric stretch of the molecule, while (Q x ,Qy) or (p,4>) 
describe bends and the asymmetric stretch. These coor- 
dinates are very useful here for our Taylor expansion of 
the body-frame quantum defects around an equilibrium 
position. We use the following forms, 

Moo (Q) = Moo (Q = 0)+a 1 Q 1 +a 2 Q 2 1 +a 3 Ql+a 4 p 2 , (31) 

Mii(Q) = M-i-l(Q) (32) 
= Mil (C = 0) + &1Q1 + b 2 Q\ + b 3 Ql + 8 P \ 

and 

M1-1 (C) = M-n (C) = Ap. (33) 

The form of our off-diagonal matrix elements (Q) 
and (Q) differs from the usual adopted form in Rcf. 
[|| by a phase factor exp (±i<fi). This different phase con- 
vention is due to the fact that the usual adopted form was 
derived in a body frame that is rotated from our body 
frame by 4>/2. We develop the detailed proof in Apendix 
B. The effective quantum numbers are calculated by di- 
agonalizing the quantum defect matrix. This calculation 
gives, 

y n ,^ (C) = n - [Mil (C) ± iMi-i (C)|] , (34) 

an expression that can be used to fit the effective quan- 
tum numbers calculated ab-initio in Ref. [20]. 

B. 3pi energy levels of H3 



TABLE II: A comparison between several of our calculated 
3pi H3 energy levels with empirically fitted experimental en- 
ergy levels 



Label" 


E ca i b 


Efu c 


Differences d 


N,g,U 


(cm-) 


(cm-) 


(cm-) 


0,1,1 


12967.8 


12966.863 


0.9 


1,0,1 


12999.1 


12998.196 


0.9 


1,1,1 


13052.3 


13050.966 


1.3 


1,2,1 


13066.9 


13068.700 


-1.8 


2,0,-1 


13139.9 


13138.608 


1.3 


2,1,-1 


13056.1 


13056.588 


-0.5 


2,1,1 


13221.1 


13219.125 


1.9 


2,2,1 


13234.3 


13235.522 


-1.2 


2,-3,1 


13203.8 


13212.055 


-8.3 


3,0,1 


13450.7 


13446.072 


4.6 


3,1,-1 


13300.6 


13300.119 


0.5 


3,2,-1 


13160.5 


13165.030 


-4.5 


3,2,1 


13485.4 


13483.545 


1.9 


3,3,1 


13453.2 


13460.934 


-7.7 



a The label denotes the values of N, G, U adopted in Ref. 0] 

to fit the experimental energy levels. 

b Theoretical results calculated in this study. 

c Empirical fits for experimental energies determined in Ref. 

0- 

d Differences theory - experiment. 



where the explicit expressions can be found in Ref. [l| 
and Ref. Q . Table 2 compares our MQDT result with 
the experimental energy levels calculated from the fitted 
effective Hamiltonians. The labels N, g, U are fitting pa- 
rameters, where N is also the total angular momentum 
of H3, and g is related to the quantum number G by 
G = \g\. Evidently our calculations are in good agree- 
ment with the fitted and recalculated experimental re- 
sults, with differences of around a few cm -1 . 



We calculate the 3pi energy levels of H3 and compare 
them with empirical fits from Ref. Q- Quantum de- 
fect parameters in Eqs. (|3"T1l3"3"l) arc extracted from the 
ab-initio calculation in Ref. |20j. To fit the experiment 
results, we shift the quantum defects at equilibrium posi- 
tions ^00 (Q = 0) = 0.0683 and p lx (Q = 0) = 0.4069 by 
a small amount, 0.0043 and 0.0021 correspondingly. 

In Ref. 0] , Vervloet and Watson studied the H3 emmi- 
sion lines of (3s, 3po, 3c?) — > 2po bands and (3s, 3po, 3d) — > 
3po bands. They then fitted the lines with effective 
Hamiltonians of the following form, 

BN (N + 1) + (C - B) K 2 
-D N N 2 (N+1) 2 
-D NK N{N + l)K 2 
-D K K 4 + ... 



IV. HIGHER ANGULAR-MOMENTUM STATES 

For higher electronic angular-momentum states with 
I > 1, the orbits are nonpenetrating and the short-range 
interaction is negligible. The long-range multipole po- 
tential model employing perturbation theory has success- 
fully described the high orbital angular Rydberg states 
of H2 [22j, |23| . In this work, we include the perturbations 
and interactions between levels of different n (principle 
quantum number) and I (angular-momentum quantum 
number) in a systematic fashion by incorporating the for- 
malism of MQDT [24[ ■ We use this long-range model to 
calculate the Rydberg states of H3 with I > 2. 

For a Rydberg electron with high orbital angular mo- 
mentum (I > 2 in the case of H^"), the effects of core 
penetration are negligible. Hence, the interaction be- 
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TABLE III: Comparison between several of our calculated 3d 
energy levels of H3 with experimentally-determined energy 
levels [J. 



Label" 


E C al b 


Efu c 


Differences d 


N+, K+,N 


(cm" 1 ) 


(cm" 1 ) 


(cm- 1 ) 


2,1,0 


17399.14 


17415.86 


-15.89 


2,2,0 


17058.41 


17039.61 


18.80 


1,0,1 


17284.81 


17296.57 


-11.76 


3,0,1 


17742.32 


17741.29 


1.03 


1,1,1 


17005.99 


16991.72 


14.28 


2,1,1 


17403.89 


17412.83 


-8.94 


3,1,1 


17698.40 


17700.43 


-2.02 


2,2,1 


17107.24 


17094.12 


13.13 


3,2,1 


17540.96 


17557.32 


-16.35 


3,3,1 


17204.46 


17188.48 


15.98 


1,0,2 


17011.36 


17001.08 


10.27 


3,0,2 


17643.36 


17655.58 


-12.21 



a The label denotes the quantum numbers N + , K + , N . 
^Theoretical results calculated in this work. 
^Empirical fits to the experimental energies determined in 
Ref. 0. 

d Differences between theoretical and experimental results. 



hence will be treated in perturbation theory. We also find 
that the quantum defect is small (of the order of 0.01), 
and the coupling between vibrational ground states of 
H^" to excited vibrational states are negligible. Hence, in 
the rovibrational transformation, we only include the vi- 
brational ground state. The body-frame reaction matrix 
thus can be written as 

K AA , n-nf drf nl (r) (Y lA \ V eff \Y W ) f nl (r), (38) 

where /„; is the regular Coulomb function with I = 1 
as the angular momentum quantum number and n as 
the principal quantum number. As the quantum defect 
for (i-wave electrons are small, we can use integers for n 
in calculating the radial functions. Here, r is the elec- 
tronic radial coordinate. Again, we perform a rovibra- 
tional transformation (with only the vibrational ground 
states) to obtain the laboratory-frame K matrix and fi- 
nally, calculate the energy levels using the standard de- 
terminantal equation of MQDT. 

Table 3 compares our theoretical calculations with 
the experimentally-determined 3d energy levels Q- The 
agreement somewhat poorer than the p-wave case. 

V. SUMMARY 



tween the Rydbcrg electron and the ion core can be ap- 
proximately described by two effects. First, the inter- 
action potential between the Rydberg electron and the 
molecular ion is expanded into a multipole series, where 
the quadrupole moment of the Hjj" core is the leading 
anisotropic term. Second, the induced dipole moment 
of the ion core interacts with the Rydberg electron by 
a potential characterized by the polarizability of the 
core. All higher angular momenta and higher-order po- 
larizabilities are neglected here, as well as the anisotropic 
portion of the polarizability interaction. 

In this approximation, the Hamiltonian is given in 
atomic units by 

H = -iv 2 -- + V ef f + H core , (36) 
2 r 

where H core is the rovibrational energy of the core. 
V e ff includes quadrupole and polarizability interactions: 

V ef f = V quad + V pol = -^P 2 (co S 0O-^-^2 (COSO • 

(37) 

where Q2, a, and 7 are respectively the quadrupole 
moment, isotropic polarizability and the cylindrically- 
symmetric anisotropic polarizability. Other components 
of the quadrupole moment tensor vanish for the undis- 
torted equilateral triangle configuration. For the vibra- 
tional ground state, Q2, a, and 7 are taken from table III 
of Ref. [25| . The polarizability and quadrupole interac- 
tions arc much smaller than the Coulomb interaction and 



In this work, we have calculated the Rydberg en- 
ergy levels of H3 molecules. Using an accurate ab-initio 
quantum-defect surface and ab-initio core energies of , 
our theoretical results for the p-wave Rydberg states 
from the present MQDT calculations are in good agree- 
ment with experimental results from J. K. G. Watson 
[?J ■ We also study higher- momentum states by a using a 
long-range multipole potential model in conjunction with 
MQDT, and find encouraging agreement with experimen- 
tal results from Ref. Q- 
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Appendix A: Permutation symmetry of the basis 
functions 

For convenience, we first consider trial basis functions 

as 

<E> try = e^TZf +m+ (a, (3, 7) ^ (9) , (Al) 

and use them to construct basis functions with proper 
permutation symmetry Eqs. (|12a| - |12b[) . The continuity 
condition for Smith- Whittcn hyperspherical coordinates 
[It], Hil requires that, 

$try {0, <p, a, p, 7) = $try (0, V + 2tt, a, /3, 7 + tt) . (A2) 
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TABLE IV: Permutation symmetry for basis functions in the 
different degrees of freedom. 



Permutation 
Operation 



P12 

P23 

P31 
P12P31 
P12P23 



K 



K+m+ 



sN' t +K' 1 



TZ 



(-) 
(-)"" 



( _ } iV+ ^iV+ 



if +771+ 



if+m+ 
JV+ 



(-)' 



K "K+m+ 



V N+ 



i27t 9j /3 ( j > 7 
i2irg I x.I 
i2 7 r 9/ /3 < j > / 

e i4,w/3gj 



where n is any integer. Under the conditions Eqs. dASal 
IA5b[) . the factors in the first line of the right hand side 
of Eq. (|A"4f vanish if m 2 = 0, gi = 0, K + = when 
is even. Therefore, the basis functions that obey the 
permutation symmetry are: 

= «i (*) ^"^U^ , (A6) 
if m 2 = 0, A" + = 0, gi = 0, and N + is odd, otherwise, 

V2 



f«j (9) 



(-f ++K+ e 



(A7) 



This boundary condition leads to the condition that 
K + /2 + ?ri2 must be integral. Hence, if A" + is even, 7712 
must be integral; if K + is odd, ttt.2 must be half integral. 
We remark here that the parity II + is determined by K + 
only; II + = +1, if K + is even, and LT+ = —1, if K + is 

odd [nm. 

The permutation symmetries for the basis functions 
chosen for each degree of freedom are shown in Table 4. 
Here, & is the nuclear-spin basis function defined as in 
Ref. [Hi- gj equals zero for the ortho state, and ±1 for 
the para state. The rotational part 7w+ + (a,/?, 7) is 
given by, 



n 



N+ 

K+m+ 



(a, &7) 



2N+ + 1 



8tt 2 



D 



N+ 

m+K+ 



(a, p, 7) 



(A3) 

where D^ +K+ are the Wigner D functions of the Eulcr 
angles. The phase of the Wigner function is chosen as in 
Varshalovich et al. [27| ■ N + is the total angular momen- 
tum of the ion, K + is the projection of N + onto the body 
frame z-axis, and m + is the projection onto the labora- 
tory frame Z-axis. We also use a set of fifth-order basis 
splines Uj (9) to expand the wave function in 8. Since 
Uj (9) is unaffected by permutations, it is not shown in 
Table 4. 

Application of the antisymmetrization operator A in 
Eq.([l3]) to Eq.flSI]) leads to: 



try 



(9) 1 + (-) A+ e i ^ L ^ m2+9/) + e'^^+si) 



5 if(mj+jr) (-) N+ 



-imif,^N+ fol 



(A4) 



Following the fact that 777,2 is integral (half integral) if 
K + is even (odd), and that gj equals or ±1, we see 
that the factors in the second line of the right hand side 
of Eq. (|A4[) vanish unless, 



m 2 + gi = 3tt. if K iseven, 



m 2 + gi = 3tt. + 3/2 if K + isodd, 



(A5a) 



(A5b) 



where 7772, gi and K + obey the relations Eqs. (|A5al 

mi). 

The explicit form of ^ j^N+K+A m ^ci- P^l i s closely 
Insertion of Eq.(|M]l into Eq.([2"T]) 



related to <f> 
yields 



N + m + gj 



jm,2N+ K+ A 



u 3 (9) (-) 



\l—A im^if 



-g I ^ J l,-A;N,K lx -Krm 



(A8) 



with K = K+ + A, if 7773 = 0, gi = 0, K+ = and N+ 
is odd. Otherwise, 



foNmgj 
^jm 2 N+K+ A 



u 3 i ) / \i— A 

e * 9 7°I,-A;IV,K'%m 



(-)' 



~9i L—A;N,K Km 



where K 



-K+ + A. 



(A9) 



Appendix B: Body-frame quantum defect matrix 
elements 



The form of reaction matrix K describing the Jahn- 
Teller coupling of the p-wave electron can be written as 



K AA < = I K 10 K n Kx-i 



(Bl) 



In perturbation theory, the matrix elements of the reac- 
tion matrix near the equilibrium position obey, 

Kaa> « -tt /" dr'f e iA (r) Uw (r) (B2) 
x J dw'X K {6' e , <p' e )AV e (Q; r e ') Y lA , (6' e , ip' e ) , 

where AV e - V e (Q; r e ') - K (Q = 0; r e '), and V e (Q, <) 
is the interaction between the ion and the Ryd- 
bcrg electron excluding the Coulomb potential, Q 
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are the vibrational symmetry coordinates, and r e ' = 
{r', 9' e , <p' e } are the spherical coordinates of the elec- 
tron in the body-frame. f £ i\ (r 1 ) is the regular phase- 
renormalized Coulomb wave function defined by, f e i\ = 
fel cos (npAA (Q = 0)) - g e i sin (npAA (Q = 0)), where 
{/ei;Sei} are the usual Coulomb wave functions with en- 
ergy e and angular momentum I. /jlaa (2 = 0) are the 
constant zero-order term of the diagonal quantum defect 
matrix elements. 

To explore the symmetry properties of Kaa' , we write 
the Taylor expansion in terms of Q = {Qi, Q+, Q-}, 
where Q+ = pe+'^ and Q- = pe~^, 

Kaa> = A AA' Q» + \Y, ^aa' Q»Qh + — ( B3 ) 

where each summation is over the subscripts 1, + and 
— . However, the expansion Eq. (|B3|) is not valid in 
our body frame coordinates because AV e (<2;r e ') is not 
single-valued at Q + = Q_ = 0. This fact is demonstrated 
in the following discussion of our body frame coordinates. 

The body frame coordinates used here are the same as 
the coordinates defined in Eq. (7) of Ref. [l7| . but we 
denote them as x'y'z' instead of xyz to make our nota- 
tions consistent. The z'-axis is perpendicular to the plane 
defined by the three nuclei and the z'-axis is associated 
with the smallest moment of inertia. After manipulating 
Eqs.(2), (3), and Eq.(7) of Ref. Qjj with some algebra, 
we can write down the cartesian coordinates of the posi- 
tions of the ith nuclei (x^, y[, z[) as 

^ = ^003 ^-1)003(1+^), (B4a) 

z[ = 0. (B4c) 

where #1 = 57r/6, $2 = — 7r/2 and #3 = tt/6. When 9 is 
very small, following Eqs. <j3j> and Eqs. ([30]) the two sets 
of coordinates {Qi, p, 4>} and {R, 9, tp} have the following 
relationship, 

Q 1 = 3 1 / 4 / (R - Ro) , (B5a) 
p = 3 1/4 fR9/2, (B5b) 

= ip - 2tt/3. (B5c) 

where Ro = 3 1 / 4 r cqui . Therefore, when Q + = Q_ = 0, 
and hence p = 0, the positions of the ith nuclei (x\, y[, z'j) 
can be written as, 

x' i (Q 1 ,Q+ = 0,Q. = 0) = 

£(* + i&7)«(! + * + §)- (Bfe » 



y' i (Q 1 ,Q+ = 0,Q-=0) = 




Zi(Qi,Q+ = Q,Q- = 0) = 0. (B6c) 

These equations show that the positions of the three 
nuclei can not be expressed by Q±, Q+ and Q- when 
p = 0, and hence the expansion Eq. (|B3|) is not valid as 
AK(Q;r e ') is not single-valued at Q + = Q_ = and 
therefore is not infinitely differentiable. However, if we 
define another set of x',y', z' axis by a passive rotation 
through </>/2 about z'-axes, the positions of nuclei have 
defined values when Q + = Q~ = in this new frame: 

x' i (Q u Q+ = 0,Q- = 0) = 

yl(Qi,Q+ = 0,Q- = 0) = 

z^(Q 1 ,Q+ = 0,Q- = 0) = 0. (B7c) 
Hence, the following expansion is valid, 

K AA > = E ^AA' Qv + \Y1 ^AA' QvQ» + ■ ■ ■ , (B8) 

and therefore using the analysis of Longuet-Higgins [28| 
in which A'aa is expanded to third order in Qi and second 
order in p, gives: 

Kaa = kIaQi + KaaQI + ^aa<9? + ^aa~ p 2 . (B9) 
We also have 

^aa±i = Aa±ia = 0, (B10) 

and 

^aa±2 = ^a±2A = K\ K±2 pe^, (Bll) 

to first order in p. The matrix elements of Kaa' also 
obey 

^aa' « -tt J dr'f £lA (r') Uw (r') (B12) 

x J d£X (e' e ,^ e )v e (Q, r e ') Yia> , 

where we use the tilde notation to stress that these are 
expressed in terms of coordinates in the new frame. The 
rotation to the new frame from our original one has the 
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following effects: 6 e = 8 e and ip e — ip e — 0/2. Hence a 
comparison with Eq. (|B2[) gives 

i?AA'=^AA'e< A - A ')*/ 2 . (B13) 
Because the reaction matrix K and quantum defect ma- 



trix fi are related by /i = fj, (Q = 0) + arctan(i4T) /i w 
A* (Q = 0) + A'/tt + O (A 3 ) /n + at least to first or- 
der, the K matrix and \i matrix have the same symmetry 
properties, and hence we write the form of /x matrix as 
Eq.® and Eas. (pT1l33j) . 
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